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O The quantum Lifshitz model provides an effective description of a quantum critical point. 

> It has been shown that even though non-Lorentz invariant, the action admits a natural su- 

k3l persymmetrization. In this note we introduce a perturbative framework and show how the 

^ supersymmetric structure can be used to greatly simplify the Feynman rules and thus the 

^ study of the model. 



Abstract 



1 Introduction 



In 1941, Lifshitz 111 introduced models with anisotropic scaling between space and time in 
the context of tri-critical models. Since then, such models have been studied in the context 
of solid state physics. Materials with strongly correlated electrons, such as copper oxides, 
show this type of critical behaviour, and also the smectic phase of liquid crystals for example 
can be described this way. Our treatment is based on quantum Lifshitz models as were 
studied in [Zf. Quantum Lifshitz points are especially interesting, since they are quantum 
critical points O, i.e. points at which a continuous phase transition happens at T = which 
is driven by zero point quantum fluctuations. A quantum Lifshitz point is characterized 
by the vanishing of the term (V^)^ in the effective Hamiltonian. While scale invariance 
is conserved, this gives rise to an anisotropy between space and time. This anisotropy is 
quantified by the dynamical critical exponent z, 

t \H, X Ax. (1.1) 

For models in 2 + 1 dimensions at a Lifshitz point, z = 2, as opposed to the Lorentz invariant 
z = 1. 

Models at a Lifshitz point have recently met with a large amount of interest beyond 
their original field of applicatiorj^ A 3 + 1 dimensional theory of gravity with z = 3 put 
forward by Hof ava [4| has generated a big echo. But also in the context of the AdS/ CFT 
correspondence, interest in gravity duals of non-Lorentz invariant models has arisen, see 
e.g. mmiZllHl- In IZl in particular, a gravity dual for a Lifshitz type model with z = 2 was 
proposed. As discussed recently in ||9|, it seems difficult to find string theory embeddings 
for gravity duals of Lifshitz-type points. 

While often, calculations are easier to do on the gravity side of the correspondence, we 
are able to perform a number of calculations directly on the field theory side, which are 
presented in this article. Apart from being of interest directly for statistical physics, our 
results can serve as a point of reference for comparison to results derived on the gravity 
side. 

In IITOl it was shown that systems of Lifshitz type in (d + 1) dimensions admit a natural 
supersymmetrization, a property which results from their relation to d-dimensional models 
via a Langevin equation. The quantum Lifshitz model in (21, described by the action 

S[<p] = I dtdx f + (did'cp)^ , x = Xi, f = 1,2, (1.2) 



^We cannot give an extensive account of the existing literature here and thus content us to only mention a 
few main examples. 
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can be thought of as descending from a free boson in two dimensions with action 



W[cp] 



dx 



(1.3) 



This formulation allows the generalization of the quantum Lifshitz model to massive and 
interacting cases. It becomes possible to consider the class of models satisfying the detailed 
balance condition whose (bosonic part of the) action takes the form 



S[<P] = J 



dtdx 



5<p 



(1.4) 



where W[^] is a local functional of the field (p{t, x). The structure due to the Langevin equa- 
tion implies supersymmetry in the time direction, so that the complete action includes also 
a fermionic field. It is given by 



S[(p,-ip,-ip] = J dtdx 



d 5'^W\<p\\ 



dt 



(1.5) 



This is the supersymmetric theory we focus on in this work. 

A major advantage of models with this structure is that they can be studied very effi- 
ciently by using a perturbative expansion of the underlying Langevin equation, as proposed 
in |[TT|. In this way, the cancellation of bosonic and fermionic terms in the perturbative ex- 
pansion becomes automatic. In consequence, there is a great simplification of the Feynman 
diagrams of the theory in + 1) dimensions which are reformulated in terms of those of the 
d-dimensional system described by W[^], plus a set of additional rules. If we consider only 
n-point functions for the bosonic field (p, all the fermionic contributions are automatically 
accounted for, so that it is not even necessary to introduce a fermionic propagator. 
For relativistic theories, this construction is possible only for d = and d = 1. Giving up 
Lorentz invariance, we concentrate on d = 2, which - as we show in the following - is the 
critical case. The generalization to any d is however clear. 

In the following we derive 

• the expression for the propagator of the free Lifshitz scalar (Sec. |3.1| ; 

• the Feynman rules for the simplest generalization to the interacting case (Sec. |3.2|; 



a scheme for UV regularization (Sec. 3.4 1 



As examples, the three-point function (Sec. 3.3.1 1 and the one-loop propagator (Sec. 3.3.2| 
are discussed. 
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2 The Langevin equation and the Nicolai map 



Having chosen to study the supersymmetric extension of the quantum Lifshitz model, we 
can make use of the Nicolai map | [T2| . In a supersymmetric field theory, a Nicolai map is a 
transformation of the bosonic fields 



[f,x) rj{t,x) , 



(2.1) 



such that the bosonic part of the Lagrangian is quadratic in tj and the Jacobian for the trans- 
formation is given by the determinant of the fermionic part: 



Sb 



dtdx 



rj{t,xf 



det 



5rj 
3^ 



VipVtp exp[-SF] 



(2.2) 
(2.3) 



Following flT], we would like to interpret the mapping in Eq. (2.1 1 as a Langevin equation 
for the field (p{t,x) with noise t]{t,x). More precisely, we want to show the equivalence of 
the action in Eq. (|1.5| to the Langevin equation 



d(p{t,x) _ 5]N 



+ f]{t,x). 



dt 5(p 

The correlations of rj, which is a white Gaussian noise (as in Eq. ( |2.2| ), are given by 

(?/(t,x)) =0, {rj{t,x))r]{t' ,x!)) =15{t-t')5{x-^) . 



(2.4) 



(2.5) 



A stochastic equation of this type, where the dissipation term depends on the gradient of a 



function of the field is said to satisfy the detailed balance condition. Equation (2.4 1 has to be 
solved given an initial condition, leading to an 7/-dependent solution (p,^{t,x). As a conse- 
quence, (p,j [t, x) becomes a stochastic variable. Its correlation functions are defined by 



{(pfj{tl,Xl) . ..(p,^{tk, Xk))^ 



JVf] exp [-If dtdxr]^{t,x)] (p,j{ti,xi) . ..(p^{ti,Xk) 
J Vt] exp [-If didxf/2(f,x)] 



(2.6) 



A necessary condition for this approach to work is that thermal equilibrium is reached for 
t ^ CO, and that 



lim {(p,j{t,xi) ...(p,^{t,Xk)) = {(p{xi) ...(p{xk)) , 



(2.7) 



i.e. that the equal time correlators for tend to the corresponding quantum Green's func- 
tions for the d-dimensional theory described by W[^]. Since (p{t,x) is a stochastic variable. 
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the expectation value of any functional F[(p] is obtained by averaging over the noise: 

(f[</>]), = ^J VfjF[<p]e-y^'^^'^^''^^\ (2.8) 
where the partition function is defined by 

2 = J P;ye-2/<ifdx,?(f,x)\ (2.9) 



It is convenient to change the integration variable from t] to (p. The expression becomes 

21 



Z = j V(p det 



exp 



dtdx 



. 5W\ 



(2.10) 



The Jacobian can be expressed by introducing two fermionic fields tp{t,x) and ^{t,x) such 



that (as in Eq. (2.3|) 



det 



VtpVxp exp 



dtdx 



3^W 



(2.11) 



In this way we can directly read off the (d + 1) -dimensional action that, up to a boundary 
term, reproduces Eq. ( |1.5| : 



S[^,xp,xp] = J dtdx 



(2.12) 



Finally, one can show [10| that in a Hamiltonian formulation, the supersymmetric ground 
state is the bosonic state 



|Yo) 



(2.13) 



as is already the case in the standard quantum Lifshitz model. 

This construction bears an obvious resemblance with the stochastic quantization of the 
d-dimensional theory described by W [(p] . We would however like to stress a fundamental 
difference. In the case of stochastic quantization, one is only interested in the f — > oo limit 



and hence in the ground state. This means that the action in Eq. ( 1.5 1 is seen as a topological 
theory. Here, on the other hand, we wish to study the finite-time behaviour of the system, 
and the same action is taken to describe a conventional supersymmetric model. 

In the following, we will work in d = 2, but in principle, all calculations are equally valid 
for general d and z = 2. The case d = lis arguably the most interesting, since it corresponds 
to the Lifshitz point with its quantum critical behaviour. 
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3 Perturbative Solution of the Langevin equation 



In the following, we will show how to perturbatively solve the Langevin equation ( |2.4| , 
which gives rise to the dynamics of the Lifshitz model. As we will see, the main advan- 
tage of this approach is that the perturbative expansion is realized in terms of the Feynman 
diagrams for the theory in d dimensions which does not include fermionic contributions. 

To solve a transport equation like the Langevin equation in Eq.( |2.4| |, it is convenient to 
consider an integral transform. The choice of transform depends on the choice of boundary 
conditions. In space, the natural choice is given by requiring the field to vanish at infinity, 

(p{t,x)- .0. (3.1) 

|x|^oo 

This means that the Fourier transform is well defined: 

(p{t,k) = fdx [e-*"(/)(f,x)l . (3.2) 



For the time direction, we have two possible choices: 

(a) The field vanishes at negative infinity. If we impose (p{t,x) > 0, we can define a 

t — > — CQ 

Fourier transform in time and use 

/CO 
dt [e-"^^(p{t,k)] . (3.3) 
-CO 

(b) Initial condition at f = 0. If we impose (p{0,x) = (po{x), it is convenient to define a 
Laplace transform in the time direction: 

/■CO 

(/)(s,k) = 1^ dt [e-''cp{t,k)] . (3.4) 

Note that the first choice preserves time-translation invariance which in case (b) is broken by 
an extra mode that describes the evolution of the initial condition. On the other hand, if the 
kernel of the Langevin equation ^ is positive definite (as it is for the cases we are considering 
here), this extra mode decays exponentially, and the large-time behaviours of both choices 
coincide. In other words, one can without loss of generality choose to impose the initial 
condition (b) and then take the large-time limit to recover the finite-time Fourier transform 
behaviour of case (a]|^ From now on, we will consider the Fourier-Laplace transform {i.e. a 
Fourier transform in space and Laplace transform in the time direction): 

<^(s,k) = j^dt j dx [e-*"-'^</)(i,x)] . (3.5) 



^An analogous problem was solved by Landau in 1 13 1 in the study of oscillations in plasma. 
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3.1 Free propagator 



As a first application, let us consider the action obtained by adding a relevant perturbation 
to the quantum Lifshitz model, described by 



S[(p,ip,ip] = dtdx 



1 



(p^ + {did'(py + m%(pd'(p + m^<p^ + fermions 



(3.6) 



According to the argument in Sec. |2| this is equivalent to the Langevin equation corre- 
sponding to the massive boson in 2 dimensions described by the functional 



dx 



(3.7) 



After the integral transform, the Langevin equation (2.4 1 takes the form 



(3.8) 



where we introduced = (k^ + m^). The Gaussian noise ?/(s, k) has the two-point func- 
tion 

2{ 2Tif3{k + k') 

s + s' 

The retarded Green's function for this problem is the solution to the equation 



(//(s,k)//(s',k')) = 



(3.9) 



sG(s,k) = -n^G{s,k) + l. 
It follows that the solution to the Langevin equation is 

(p{s,k) = G(s,k)77(s,k) + G(s, k)</)o(k). 



(3.10) 



(3.11) 



and since the Laplace transform exchanges point-wise products and convolution products, 
one finds that the field ^ as a function of time can be written as 

(p{t,k) = G{t,k)*?]{t,k) + (Po{k)G{t,k) = f dT [G{t - T,k)rj{T,k)]+ (pQ{k)G{t,k) . 

Jo 

(3.12) 

More explicitly, using the fact that 



G{t,k) = e 



(3.13) 



we find the solution 



(P{t,k) = dT 



//(T,k) +(po{k)e 



(3.14) 
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Having expressed (p as a function of the noise, we are now in a position to evaluate the 
two-point function: 



D(s,k;sM<') = {cp{s,k)<p{s\k')) = G(s,k)G(s',k') (j/(s,k)//(s',k')) 



2(27rr^(k + kO 
(s + Q2) (s' + Q2) (s + s') ■ ^ ' 



Taking the (bidimensional) inverse Laplace transform, this becomes 
D(f,k;f',k') = / ^ ^ e^^+^'''D(s,k;s',k') 



= (2nf'—^ _^(k + k'). (3.16) 



Two limits are interesting: 



• For t, t' — > 00, the second exponential is suppressed and the two-point function be- 
comes 

D{t,h;t',\^) > (2nf—^s{h + h'). (3.17) 

Its Fourier transform is given by 

D(..,k;.-,k-) = -"" + fS+'-'' , (3.18, 

CO +12 

as found in pj. 

• For t = t' ^ oo, the two-point function reproduces the usual bosonic propagator in 
d dimensions (as expected from the general structure of the Langevin equation and 
shown in Eq. ( |2.7[ |): 

D{t,k;t,k') > (271)2 

Note that for m = 0, this means that for large times, the propagator will behave poly- 
nomially and the theory is critical, just like in the case of a two-dimensional bosonic 
field. 



3.2 Interacting theory 



As an example of an interacting theory, let us consider the action descending from the theory 
of a massive boson with a cp^ interaction. This is described by a Langevin equation with 
functional 



w[<p] 



dx 



(3.20) 
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The action in (2 + 1) dimensions in Eq. (1.5 1 is immediately found to be 



S[(p,xp,ip] = dtdx 



1 



+ {did'cpy + m^di(p d'(p + g'^(p di(p d'(p + m^(p^ + g^m^(p^ + g'^^^+ 



+ fermions 



, (3.21) 



and in particular, for the critical m = case, 

"1 . 



dtdx 



(p^ + {did'(pf + g^(p di(p d'(p + g'^cp'' + fermions 



(3.22) 



Note that the coefficients of the three- and four-point interactions are not independent, but 
fixed by the detailed balance condition. It is in fact the detailed balance, which here plays 
the the role of a symmetry the theory has to fulfill, which keeps terms other than those given 
in Eq. ( 3.21| from appearing. This relation between the different coupling constants is a 
property which is accessible to experimental checks in materials which are described by a 
Lifshitz point effective action and thus is a testable prediction of this framework. 

The most effective way of making use of this symmetry consists once more in starting 
from the corresponding Langevin equation, which takes the simple form 



dt(pit,x) 



-V^(p{t, x) + m^(p{t, x) + g^<p^{t,x)) + ^{t,x) , 



(3.23) 



where y]{t,x) is the same Gaussian noise as in the previous example. Taking the Fourier- 
Laplace transform as above, one has to deal with the quadratic term that can be recast in the 
following form: 

-g^ j"" dt j dxe-'^ ''-'^(j){t,xf = 



g^ / £f[ki,k2,Si,S2] 



^(k - ki - k2) 



S — Si — S2 

The Langevin equation becomes an integral equation: 

"^(k-ki-k2) 



Si,ki)^(S2,k2) 



. (3.24) 



(s,k) = G(s,k)//(s,k) -^^G(s,k) / d[ki,k2,Si,S2 



^(Si,ki)(/>(s2,k2) 



S — Si — S2 

+ G(s,k)(/>o(k). (3.25) 



This type of equation can be solved perturbatively in g^, using the usual Feynman diagram 
techniques by denoting G (s, k) with an arrow and the noise // (s, k) with a cross. In particular. 
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element 



Fourier-Laplace 



Fourier 



-X- 



k k' 

S2 k2 

Si 



s-l 



S3 ks 



G(s,k) 



1 



S + Q2 



2(27r)^^(k + kO 
(s + Q2) (s' + Q2) (s + s') 

2<5(ki - k2 - ks) 

Si - S2 - S3 



G(a;, k) 



(271)^+^ ^(k + kO^(a; + a;0 
a;2 + Q4 



^^^(ki - k2 - k3)^5(a;i -CO2- co^) 



Table 1: Feynmarv rules for the cubic theory 



the field ^ can be expanded as 




+ 




+ 



(3.26) 



The Feynman rules for this model are summarized in Table [Tj Note that even if the action 
in Eq. ( |3.21 1 has two cubic and a quartic interaction, the Feynman diagrams obtained from 
the Langevin equation only have one cubic vertex. This might at first seem surprising, but 
is again a simplifying consequence of the detailed balance condition. 



3.3 Examples 

3.3.1 Three-point function 

As a first example, let us consider the three-point function (^(si,ki)^(s2, k2)(^(s3, k3)) at 
tree level. Expanding the field as in Eq. ( |3.26| , we find that at tree level the function is given 
by the sum of three contributions: 

S2 k2 S2 k2 k2 

Si 



{(p{Si, ki)</)(s2, k2)(/'(S3, ks)) 



ki 




(3.27) 



ss'ks sa'ka 

Using the Feynman rules in Table [ij it is immediate to find that each diagram gives a contri- 
bution of 



dsadst 
47r2 



(3.28) 
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where {i,j,k} take the values {1,2,3} and their cyclic permutations. After an inverse 
Laplace transform in the time direction and taking the large time limit (which together is 
equivalent to taking the Fourier transform), the three-point function becomes: 



~^ n2n2n2((n2-n2)''-nf) ' ' ' 

where tjj = tj — tj and ti < t2 < t^. A consistency check for this expression is obtained by 
considering the equal time case fi = ^2 = ^3 ^ which reproduces, as expected, the usual 



bosonic result (as in Eq. (2.7|): 



{cp{tMm,k2)cp{t,k3)) ^ ^2^2^2 ■ (3-30) 



3.3.2 One-loop propagator 

As a next example, let us consider the one-loop correction to the two-point function (^(si, ki)(^(s2, k2)) . 
Using the expansion in Eq. ( 3.26| l, we find that 



(a) k2 k2 k2 

(b) (c) (d) (3.31) 

The first term is just the usual propagator D(si, ki; S2, k2). The contributions from the dia- 
grams (b), (c), (d) are as follows: 

4 G(si, k)D{sa, ki; s^, -ki)D(sb, k2; Sc, -k2)G(s2, k) 
^ (Si - S„ - Sfc) (S2 -s^-Sc) 

4 D(si,k;Sfl,-k)D(sb,ki;Srf,-ki)G(sc,k2)G(s2,k) 

[Sc - Sa - Sb) {S2 - Sc - Sd) 

4 G(si, k)D(sa, ki; s^, -ki)G(sfc, k2)D(srf, k; S2; -k) 

(Si - S„ - Sfc) (Sb - Sc - Srf) 

The two-point function is obtained by summing up all the contributions and integrating 
over the internal momenta. Once more, the result is more transparent if we take the inverse 
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Laplace transform and consider the large time limit fi, f2 



dkidk2 



(2nt 



, , , , \ — 5Ck - ki - k2) 



, (3.35) 



where t2i = | ^2 — | • Taking the ^21 limit, one reproduces the usual bosonic two-point 
function at one loop (as in Eq. ( |2.7| ): 

/^^i 1 1 ^\ 1 , 4 /■ dkidk2 ^5(k - ki - k2) „ 

m,k)<p{t,k)) J (,^)4 040202 • (3-36) 

3.4 UV Regularization 

The two-point function we derived above suffers from an ultraviolet divergence. In order to 



regularize it, one can either change the Langevin equation (Eq. (2.4 1), or the noise correlation 
function in Eq. ( |2.5| |. Here, we follow the latter approach and show how smearing out the 
delta function in time by introducing a cut-off A results in an ultraviolet regularization in 
space which, in the large time limit, reproduces the usual Pauli-Villars result. 
Consider the noise function ?/a(^/ x) with the following two-point function: 

(//A(i,x)/;A(i',x')) =^(x-x')AV^'I'-*'I. (3.37) 

For A ^ 00, it converges to the two-point function of the usual noise: 

(//A(f,x)//A(i',x')) >23{x-x')S{t-t'), (3.38) 

rjAihx) >r]{t,x). (3.39) 

A^oo 

Applying the Fourier-Laplace transform, we get 

/ / N . / (27r)^^(k + k') a2(s + s' + 2A2) 

(,a(^,x),a(^',x')) = ^1,. (s + A2)(s- + A2) - (3-40) 

In this scheme, the Langevin equation remains unchanged, which implies that the retarded 
Green's function remains the same, 

G(.,k)=^. (3.41) 
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while the field receives a correction. In particular, for the free case: 



(pA{s,k) = G{s,k)t]A{s,k) + G{s,k)(Po{k) . (3.42) 
We are now in a position to calculate the correction to the propagator in Eq. ( 3.15| |: 



D k-.' k') - let, kU k')) - (2/rf'^(k + kO A^(s + s- + 2A^) 
DA{s,k,s,k)- {cpA{s,k)cpA{s ,k )) - ^ ^2) (g, + 02) + g,) + (s' + A2) • 

(3.43) 

To see how this corresponds to an ultraviolet regularization, let us perform the inverse 
Laplace transform: 



Q2 (A4 - Q4) 



For large times, only the first two terms in the sum contribute, and for t = t' oo, we find 
the usual Pauli-Villars propagator for the 2-dimensional boson: 

D((.k;(,k') ^ (2.)^^(k + k-) „^/^„^ . (3.45) 
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